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We propose a new method to discover light top squarks (stops) in the co-annihilation region at 
the Large Hadron Collider (LHC). The bino-like neutralino is the lightest supersymmetric particle 
(LSP) and the lighter stop is the next-to-LSP. Such scenarios can be consistent with electroweak 
baryogenesis and also with dark matter constraints. We consider the production of two stops in as- 
sociation with two 6— quarks, including pure QCD as well as mixed electroweak-QCD contributions. 
The stops decay into a charm quark and the LSP. For a higgsino-like light chargino the electroweak 
contributions can exceed the pure QCD prediction. We show the size of the electroweak contribu- 
tions as a function of the stop mass and present the LHC discovery reach in the stop-neutralino 
mass plane. 



I. INTRODUCTION 

The Large Hadron Collider (LHC) is currently col- 
lecting data at y/s = 7TcV and it is assumed that 
the integrated luminosity will reach lfb -1 next year. 
If low-energy supersymmetry (SUSY) [l[ is realized, 
detection of light supersymmetric particles may be 
around the corner 

The scalar top (stop) within the minimal supersym- 
metric standard model (MSSM) [l[ is naturally one of 
the lightest SUSY particles (sparticles). On the one 
hand, the large top Yukawa coupling leads to large 
mixing between left- and right-handed stops, reduc- 
ing the mass of the lightest stop mass eigenstatc. On 
the other hand, the large top Yukawa coupling reduces 
the stop mass at the electroweak (EW) scale via renor- 
malization group equation (RGE) running Q. The 
lightest stop, ?i, is mostly right-handed, since is 
not increased by SU(2) gaugino loops, and is more 
strongly reduced by the top Yukawa interaction Q. 

A light stop with a mass of mi < 125 GeV is vital 
for successful EW baryogenesis within the MSSM 0- 
It allows for a strong first order phase transition, 
which prevents the generated baryon asymmetry of the 
universe flfjl. [Tllj from being washed out. Furthermore, 
CP violation is needed in order to generate the baryon 
asymmetry [HJ. The Standard Model (SM) CP violat- 
ing Dirac phase is too small whereas the MSSM 
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contains additional CP-violating phases sufficient for 
EW baryogenesis 0j- 

However, new CP violating phases are strongly con- 
strained by the non-observation of electric dipole mo- 
ments ruling out large regions of the MSSM consistent 
with EW baryogenesis p, HI, ||| [l4| . Of the remaining 
MSSM parameter space, scenarios with bino-driven 
baryogenesis are probably the most promising ones 
@, Here, CP violation in the bino-higgsino sec- 
tor accounts for successful EW baryogenesis. In this 
scenario the (stable) lightest supersymmetric particle 
(LSP) is a bino-like neutralino and the Higgs mixing 
parameter fi needs to be of the order of the bino mass 
Mi. In addition, all sfermions (beside the light stop) 
are quite heavy, in order to suppress electric dipole mo- 
ments and to fulfill present bounds on the Higgs mass 
[3, E3|- The light stop should be predominantly an 
SU(2) singlet (i.e. "right-handed") in order to sup- 
press stop loop contributions to the electroweak rho 
parameter [l5| . As we show below, many of these sce- 
narios fall into the parameter space that we investigate 
in this work (49j . 

In the MSSM, the lightest neutralino, Xi> is a 
promising dark matter (DM) candidate if it is the LSP 
[la ]. However, large regions of the MSSM parameter 
space are disfavored due to a too large relic density 
of the X? Q3 • The relic density is determined by the 
thermally averaged cross section, which includes anni- 
hilation and co-annihilation processes. Other sparti- 
cles with masses not far above that of the LSP can co- 
annihilate with the neutralino and/or enhance t— or 
u— channel exchange contributions to the annihilation 
processes reducing the DM density in the universe to 
a level consistent with cosmological observations [ll| . 
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If the mass splitting between the neutralino and 
another particle is < 20%, the co-annihilation dia- 
grams are significant We consider such scenarios. 
We assume a relatively light neutralino LSP, a stop 
next-to-lightest supcrsymmetric particle (NLSP) and 
a light higgsino-like chargino xt- AH other sparticles 
are assumed to be heavy. We take neutralino-stop and 
stop-chargino mass differences of a few tens of GeV. 
These scenarios can be consistent with DM constraints 
[H, [l9[ , i. e. the thermal Xi renc density is equal to or 
lies below the observed DM density. 

For most of our scenarios the (thermal) Xi re li c 
abundance lies below the observed one if « \M X \ 
[50j . This is because a non-negligible higgsino com- 
ponent of the Xi leads to efficient Xi annihilation into 
W + W~ and Z°Z° pairs. However, there are several 
options to make scenarios with low thermal xi density 
phenomenologically viable. The simplest option is to 
have an additional DM component, see e.g. Refs. [2lj . 
Furthermore, a non-standard cosmological history can 
lead to the right Xi density [22[ . Finally, the Xi abun- 
dance can arise from the decays of mctastablc species if 
their lifetime is so large that they decay out of thermal 
equilibrium [23j . In contrast, dilution of a (too large) 
thermal Xi relic density is more difficult to achieve (8| . 

Stop pair production might be difficult to detect 
in the co-annihilation region at hadron colliders fljjj ]. 
The decays t\ — > bx\W, with b (W) the bottom quark 
(W boson), as well as ti — > tvfb, v^ib are kincmatically 
closed. The loop induced two-body decay t\ — > cx\ 
then competes with tree-level four-body decays like 
ti — > lvibx\- The latter are strongly phase space 
suppressed for small stop neutralino mass splitting, 
so that the l oop induced decay becomes the dominant 
decay mode (24. [2o|. However, the charm quark, c, is 
soft, so that the collider signature is given by two soft 
charm jets and missing energy. 

Tevatron searches for light stops, decaying to charm 
and Xi i require a minimum mass gap between the stop 
and X? of at least 40 GeV [H, |27|. Otherwise the 
charm jets are too soft to be seen above the SM back- 
grounds. As a consequence, the Tevatron is not very 
sensitive to the coannihilation region. At the LHC, 
the detection of stop pair production is expected to 
be even more difficult due to a large hadronic activity. 
However, gluinos can decay into a stop and a top. Due 
to the Majorana character of gluinos, same-sign tops 
can be produced. The resulting collider signature can 
be probed at the LHC as long as the gluino mass does 
not exceed 900 GeV [H; see also Refs. [H. Ref. ^ 
therefore investigated QCD stop pair production in 
association with a very energetic photon or jet. The 
resulting signature is one hard photon or jet recoiling 
against large amounts of missing energy. Their results 
are independent of the gluino mass and the discovery 
reach covers stop masses up to roughly 240 GeV for an 
integrated luminosity of 100fb _1 at y 7 ! = 14TeV. 

Here we instead consider the production of two stops 
in association with two b (anti)quarks. We include 



pure 0(ag) QCD diagrams as well as leading order 
0{a 2 s aw) mixed QCD-EW contributions. The lat- 
ter are due to diagrams with an on-shell higgsino- 
like chargino and substantially increase the total cross 
section. Information about the magnitude of the EW 
contributions can in principle be obtained by measur- 
ing our process and stop pair production in association 
with a hard jet [3(j. The latter process can be used to 
determine the stop mass, which in turn uniquely fixes 
the pure QCD contribution to our process. Subtract- 
ing this from the measured cross section would yield a 
determination of the magnitude of the EW diagrams. 

These diagrams are sensitive to the higgsino cou- 
pling t\ — X\ ~ b. Therefore, the measurement of 
our process allows a test of the respective SUSY cou- 
pling relation and thus a test of supersymmetry itself 
[3ll |32| . Only a few of such tests have been proposed 
so far for the LHC [3ll - l33j and none of them addresses 
interactions which originate from the supcrpotcntial. 

The remainder of this article is organized as follows. 
In Sect. |H] we describe our process and present the 
dominant QCD and EW contributions. In Sect. lIIIi we 
first discuss the dominant background processes and 
then basic cuts for a benchmark scenario before pre- 
senting our numerical results. We show the discovery 
reach in the neutralino stop mass plane. In Sect. IIV1 
we discuss possibilities to further optimize our cuts, 
which sensitively depend on the stop neutralino mass 
splitting. We conclude in Sect. [Vj 

II. STOP PAIR PRODUCTION IN 
ASSOCIATION WITH TWO B-JETS 

We consider stop pair production in association with 
two b— jets in proton proton collisions, 

pp^htlbb. (1) 

We are interested in scenarios where the mass differ- 
ence between the lightest stop and neutralino is not 
larger than a few tens of GeV as discussed in the in- 
troduction. In such scenarios, the stop decays into the 
lightest neutralino and a soft charm jet, 

h -+ xic (2) 

Due to the small t\ — x\ mass splitting, much of the 
time the c— jets will be too soft to be useful (28|. Our 
hadron collider signature is therefore large missing en- 
ergy and two b— flavored jets. We require both hard 
jets to be tagged as 6— jets, since this greatly sup- 
presses SM backgrounds. The dominant QCD contri- 
butions are generated via gg — > iii* processes, where 
the bb pair comes from additional gluon radiation split- 
ting into bb. One of the dominant diagrams is shown in 
Fig.[TJ QCD contributions with quarks and antiquarks 
in the initial state arc subdominant. For example, for 
mj-j = 400 GeV, quark-antiquark diagrams contribute 
only about 6% to the total cross section due to the 
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small qq flux at the relevant Bjorken— x (34|; this is 
smaller than the uncertainty of our leading order anal- 
ysis. For smaller stop masses these contributions are 
even less important. Therefore, we only consider gluon 
fusion diagrams; cf. Figs. Q] and [2] 




FIG. 1: Example diagram for QCD stop pair production 
in association with two b— jets via gluon fusion. 

We also include the leading EW contributions to our 
process; one of the corresponding Feynman diagrams 
is shown in Fig. [2j Diagrams with electroweak gauge 
bosons and Higgs bosons exchange are subdominant 
and are not taken into account in our analysis. Again, 
contributions with a quark and an antiquark in the 
initial state are suppressed and we do not consider 
them. If an on-shell decay of a chargino into a stop 
and a b— quark is kinematically possible, the diagram 
shown in Figure [2] is effectively only a 2 — > 3 process, 
because the (on-shell) chargino will decay into t\b\ this 
decay will almost always be allowed if i\ is the NLSP. 
If the chargino mass is not much above the stop mass, 
this process is therefore less phase-space suppressed 
than the 2 -> 4 QCD (and EW) contributions. We 
found that for a chargino with Am = m -± — m; = 
20 GeV, EW contributions are comparable to the QCD 
contributions. 




b 



FIG. 2: Example diagram for EW stop pair production in 
association with two b— jets via gluon fusion. The chargino, 
Xi , might be on-shell. 



As motivated in the introduction we assume that 
t\ is dominantly right-handed. Hence, EW contribu- 
tions are maximal if the light chargino is higgsino-like. 
If it is wino-likc EW contributions will be suppressed 
by small mixing angles, because the wino does not di- 
rectly couple to the right-handed stop. Therefore, our 
results will not change significantly if one adds light 
winos to our scenarios as long as two-body decays of 
the stop and the higgsino into winos are kinematically 
forbidden. 

As can be seen from Fig. [51 we are sensitive to the 
Xi — b — t± coupling. We will present prospects of 
determining the xt — b — t\ coupling in a later pub- 
lication. If a signal is observed, one should also be 
able to obtain information about the mass spectrum 
of the chargino and stop sector. However, here we are 
primarily interested in the discovery reach of this new 
t\ search channel. 



III. LHC ANALYSIS 

In this section we first review the dominant SM 
background processes. Next we choose a benchmark 
scenario in the co-annihilation region, which is com- 
patible with electroweak baryogenesis (if one adds an 
additional CP-phase) . We then present kinematic dis- 
tributions and discuss our basic cuts and compare the 
size of the EW contributions to the pure QCD predic- 
tion. Finally, we show the discovery reach at the LHC 
in the neutralino stop mass plane. 



A. Backgrounds 

We only consider SM backgrounds, since we assume 
that all other colored sparticles are quite heavy so that 
SUSY backgrounds are negligible. We look for SM 
processes which lead to two b— jets and large missing 
energy. The dominant backgrounds arc 

• ti production (including all top decay chan- 
nels). Top decays will nearly always produce 
two 6— jets. Since we require large missing Et, 
at least one of the W bosons produced in top 
decay will have to decay leptonically. Note that 
this also gives rise to a charged lepton (e, fx or 
t), whereas the signal does not contain isolated 
charged leptons. 

• Z(— > vv) + bb production, i.e. Z boson produc- 
tion in association with two b— jets. The Z boson 
decays into a pair of neutrinos. If the charm jets 
in the signal are very soft, this background looks 
very similar to our signal. Fortunately it can be 
directly extracted from data. One can measure 
Z(—> e+ e~ I fi + + bb, where the Z decays into 
either a pair of electrons or muons. From the 
known Z branching ratios (BRs) one can then 
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process 


Wbb 


single top 


Zbb 


tt 


a [pb] 


84 


170 


174 


800 



TABLE I: Total hadronic cross sections in pb for the main 
SM backgrounds at yfs — 14 TeV. The cross sections were 
calculated with Madgraph apart from ti production, which 
is calculated in Ref. 13811. 



with Fast Jet2 .4.1 [43| via the kt clustering algorithm 
with R = 1.0. Our event samples are then analyzed 
with HepMC2 . 04 . 02 0] and ROOT [45[. For the sake of 
simplicity, wc keep the b— and c— flavored hadrons sta- 
ble. A jet is identified as a b— jet, if a stable b— hadron 
is found in the reconstructed jet. If not otherwise men- 
tioned, we assume a b— tagging efficiency of 60%. 



obtain an estimate for the background cross sec- 
tion. However, this procedure will increase the 
statistical error due to a smaller BR of the Z 
to charged leptons compared to the decay into 
neutrinos [35j . 

• W(— >• Iv) + bb production, where the W decays 
lcptonically. Again, this background will con- 
tain a charged lepton, and will thus resemble the 
signal only if the charged lepton is not identi- 
fied. This can happen when the charged lepton 
emerges too close to the beam pipe or close to a 
jet; moreover, identification of hadronically de- 
caying r leptons is not easy. 

• Single top production in association with a 
b— quark, e.g. gu — > tbd. The second b— jet stems 
from top decay, and the missing Et comes from 
the leptonic decay of the W boson. 

We neglect QCD dijet and trijet production in our 
analysis, since a large $ T cut should remove those 
backgrounds [32l. l36l l37j . 

Estimates for the total hadronic cross sections for 
these SM backgrounds are given in Table HI The cross 
section for the ti background has been taken from 
Ref. [Hj], which includes NLO corrections as well as 
resummation of next-to-leading threshold logarithms. 
All other backgrounds, as well as the signal, have been 
calculated to leading order using Madgraph4 . 4 . 32 

M- 

We have generated 5 x 10 6 tt, single top and Wbb 
background events, respectively, as well as 3 x 10 6 Zbb 
events. 



B. Numerical tools 

The masses, couplings and branching ratios of the 
relevant sparticles are calculated with SPheno2 .2.3 
[io| . We use the CTEQ6L1 parton distribution func- 
tions and the one-loop expression for the strong gauge 
coupling with five active flavors with Aqcd = 165 MeV 
[4lT |. Our parton-level signal and background pro- 
cesses apart from ti production are generated with 
Madgraph4 .4.32 _39j . Parton-level events are then 
interfaced with Herwig++2 .4.2 [4^] for the hadron- 
levcl simulation. We also generate our ti events with 
Herwig++, fixing the normalization as in Ref. (3^. Wc 
do not consider detector effects. Jets are reconstructed 



C. Benchmark Scenario 

In order to develop a set of cuts, we introduce a 
benchmark scenario. We work in the framework of the 
general MSSM with a light (dominantly right-handed) 
stop with TOjj = 120 GeV, which is compatible with 
electroweak baryogencsis. Our benchmark point is 
also consistent with DM constraints. The lightest neu- 
tralino is bino-like with 

m _o = m h - 20 GeV. (3) 

The lightest chargino is higgsino-like with 

m-± =mx +20 GeV. (4) 

All other sparticles are decoupled. The cross sec- 
tion for our benchmark point is given in the first line 
of Table ITT1 For comparison, we separately present the 
total hadronic cross section for the QCD+EW contri- 
butions (second column) and the pure QCD contribu- 
tion (third column) as well as the ratio of the respec- 
tive cross sections (fourth column). We also display 
the cross sections for heavier stops (first column) as- 
suming the mass relation of Eq. Q . 

We can see that the cross section decreases quickly 
with increasing sparticle masses as expected. For 
example, increasing the stop mass from 120 GeV to 
200 GeV decreases the total cross sections by roughly 
a factor often. However, as we will show in Sect. IIII El 
the significance with respect to the SM backgrounds 
will decrease less rapidly. This is because the final 
state particles will have on average larger momenta 
and therefore the events will more easily pass our cuts. 
We also observe that the EW contributions are signifi- 
cant and even give the leading contribution to the cross 
section. For the case at hand, they enhance the total 
hadronic cross section by roughly 150% compared to 
the pure QCD contribution. 

Eq. ((4]) implies that the decay ti — > x^b is not al- 
lowed. Similarly, Eq. ([3]) implies that t\ — > t\i and 
ti bj(iW decays are forbidden. As mentioned in the 
introduction, our stop then decays via a flavor chang- 
ing neutral current (FCNC) decay, BlZ(ti — > c\\) ~ 1. 
This requires that the physical t\ has a nonvanishing 
c component. Even if squark flavor mixing is assumed 
to be absent at some (high) input scale, it will be in- 
duced by electroweak one-loop diagrams, due to the 
fact that quark generations do mix. This one-loop 
process is enhanced by large logarithms [24| , i.e. it can 
be understood as describing the running off-diagonal 
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14U 
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160 


(3.0 


2.2 


2.0 


180 


3.2 


1.3 


2.5 


200 


2.0 


0.81 


2.5 


220 


1.2 


0.51 


2.4 


240 


0.83 


0.34 


2.4 


260 


0.56 


0.23 


2.4 


280 


0.38 


0.16 


2.4 


300 


0.27 


0.11 


2.5 


320 


0.19 


0.081 


2.3 



TABLE II: Total hadronic signal cross sections in pb from 
the pure QCD diagrams (third column) and from the 
QCD+EW contributions (second column) as a function of 
the stop mass (first column). The mass relation of Eq. Q 
has been assumed to hold. We also show in the fourth 
column the ratio of the QCD and QCD+EW cross sec- 
tions. All cross sections were calculated with Madgraph for 
v/i = 14 TeV. See Sect. HJ for further details. 



tifih mass. Of course, it is also possible that squark 
mass matrices are not exactly flavor diagonal at any 
scale. The stop generally decays promptly, i.e. its 
flight path is much too short to be seen experimentally. 
However, depending on the size of the t\c\\ coupling, 
the lifetime of the stop can exceed 1/Aqcd ~ 10~ s 
in which case the stop hadronizes before it decays (4(| . 
We assume that BlZ(xf — > i\b) = 1; this almost 
follows from our assumption that both (third gener- 
ation) charged and neutral sleptons are heavier than 
mi - m b . 



D. Distributions 

We present in this section kinematic distributions 
at the LHC of the background and signal for our 
benchmark scenario. We show cumulative distribu- 
tions, i.e. the expected signal and the different back- 
ground contributions are stacked on top of each other. 
Note that we show the number of events on a loga- 
rithmic scale. All distributions are scaled to an inte- 
grated luminosity of lfb^ 1 at ^/s = 14 TeV. In all 
cases we require at least two b— jets with a rapidity of 
\r]b 12 \ < 2.5. We also require the jets to have trans- 
verse momenta px > 20 GeV. We assume in this sub- 
section a b— tagging efficiency of one, but we will as- 
sume later (in Sect. IIII E[) a more realistic efficiency 
of 60% when we derive the discovery potential. No 
further cuts are applied. 

In Fig. [3j we present the number of isolated charged 
lcptons (electrons, muons) for signal and background. 




We only consider isolated leptons with p T 



lepton 



> 



number of isolated leptons 

FIG. 3: Number of isolated leptons for the signal and 
SM backgrounds assuming an integrated luminosity of 
lfb -1 at = 14 TeV. Fo r the signal we assumed the 
benchmark scenario of Sect. ITTTCl i.e. m^a = 100 GeV, 
m f = 120 GeV and m~+ = 140 GeV. The distributions 
are stacked on top of each other. 



5 GeV. A lepton is isolated, if less than 10 GeV of en- 
ergy are deposited in a cone of AR = v/A0 2 + A?y 2 < 
0.2 around the direction of the lepton (not counting 
the energy of the lepton itself). The signal process 
possesses nearly no isolated leptons, because no lep- 
tons can arise at parton level. We thus employ a veto 
on isolated leptons in Sect. IIII El This cut will effec- 
tively reduce the SM backgrounds involving leptoni- 
cally decaying W bosons. Recall that we have only 
considered the Z — > vv channel for the Zbb back- 
ground. Since our fe— hadrons are stable, we have no 
isolated leptons for the Zbb background and at most 
two isolated leptons for the top backgrounds. Later we 
will impose quite stiff cuts on the transverse momenta 
of both b— jets; leptons originating from semileptonic 
b decays would therefore not be isolated. Similarly, 
the leptons resulting from semileptonic charm decays 
will be either very soft or not isolated. Recall also 
that we require the W boson in the Wbb background 
to decay lcptonically, while we do not demand spe- 
cific decay modes for the top quarks in the single top 
and tt backgrounds. However, these backgrounds can 
only produce large missing px if at least one W boson 
decays leptonically. In that case the result for the sin- 
gle top background would be very similar to that of 
the Wbb background, while the distribution for the tt 
background would peak at ni = 1. 

Fig. @] shows the pr distribution of the hardest 
b— jet. The tt and single top backgrounds give the 
hardest b— jets, because here at least one 6— quark 
arises from the decay of a heavy top quark. The 
other backgrounds have softer b— jets, which originate 
mainly from a bb pair generated via gluon splitting. 
For the signal, the dominant QCD contribution is tit* 
production with an additional bb pair from gluon split- 
ting, whereas the main EW contributions are from 
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FIG. 4: Same as Fig. [3] but now for the pr distribution of 
the hardest fe-jet. 



FIG. 5: Same as Fig.[3j but now for the pr distribution of 
the second hardest &-jet. 



tiXib production, where the second &-quark comes 
from a resonant xt decay; cf. Sect. [TXJ 

Note that for our benchmark scenario, in the mixed 
QCD+EW contribution the b— quark from the xt de- 
cay is usually much softer than the other b— quark. 
Moreover, g — > bb splitting prefers asymmetric con- 
figurations, with one b (anti)quark being significantly 
more energetic than the other (HT| ■ As a result the 
Pt distribution of the second hardest b— jet shown in 
Fig. [5] is much softer for the signal, peaking around 
40 GeV, whereas the pr distribution of the hardest 
signal b— jet peaks around 70 GeV; see Fig. 2] The 
distribution of the second b— jet in the single top back- 
ground is also very soft, since it originates from g — > bb 
splitting in the initial state, and that for the Wbb back- 
ground is soft since here the b— jets originate from 
g — > bb splitting in the final state. In contrast, the 
main contribution to the Zbb background can be un- 
derstood as gg — > bb production where the Z boson is 
emitted off the quark line; here, and in the ti back- 
ground, the difference between the pt distributions of 
the two 6— jets is therefore smaller than for the other 
processes. 

The reason for the soft px spectrum of the sec- 
ond signal b— jet is the small mass splitting of 20 GeV 
between the t\ and the xf- I n the rest frame of 
the decaying chargino, the b— quark has 3-momentum 
= 10 GeV. If we increase the mass gap be- 
tween stop and chargino, the prfa) distribution will 
be much harder, e.g. for m~± = 210 GeV we ob- 
tain \p*\ = 70 GeV. However, increasing m-± also 
decreases the signal cross section significantly. For ex- 
ample, increasing the stop chargino mass difference 
from 20 GeV to 80 GeV reduces the total hadronic 
cross section from 19 pb (see Table [TTJ) to 10 pb. In 
this case 75% of the cross section would come from 
pure QCD contributions. This reduction of the cross 
section ovcrcompensatcs the gain of efficiency due to 
the harder prfo) spectrum, i.e. the signal cross sec- 
tion after cuts also decreases with increasing chargino 
mass, although less quickly than the total cross section 




50 100 150 200 250 300 350 400 



ptmiss [GeV] 

FIG. 6: Same as Fig. [3] but now for the missing transverse 
momentum distribution. 



before cuts does. 

With the help of Figs. [5] and [SJ we found lower cuts 
of 150 GeV (50 GeV) on the hardest (second hardest) 
b— jet helpful in order to increase the signal to back- 
ground ratio. We will employ these cuts in the next 
subsection. A similar cut on the hardest jet (and on 
the missing energy) is also required by the trigger [47j . 

The missing transverse momentum, p T distributions 
of the signal and the backgrounds are given in Fig. |B1 
We observe that the signal fl distribution falls off less 
quickly than the background distributions. Therefore, 
we will employ a lower cut on p of 200 GeV. On the 
one hand, this cut increases the significance of the sig- 
nal over the top, Zbb and Wbb backgrounds. On the 
other hand, we expect that it suppresses pure QCD 
backgrounds like dijet and trijet production to a neg- 
ligible level [HHl. 

In addition to kinematic distributions we can also 
employ the number of charged particles (mainly 
hadrons) to distinguish signal from background events. 
The respective distributions are given in Fig. where 
only charged particles with pr > 2 GeV and \t]\ < 2.5 
are included. Moreover, we have assumed the W bo- 
son in single top events as well as at least one W boson 
in ti events to decay leptonically in order to obtain a 
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multiplicity 

FIG. 7: Same as Fig. 02 but now for the multiplicity of 
charged particles with \rj\ < 2.5 and pr > 2GeV; more- 
over, we have required that single top and ti events contain 
at least one leptonically decaying W boson. 

significant amount of missing energy [52j ■ The num- 
ber of charged particles, A^ c harged) m tt events is nev- 
ertheless on average larger than for the signal. The 
second W boson will usually decay hadronically, pro- 
ducing jets that are usually considerably harder than 
the c— jets in the signal. Recall that single top pro- 
duction is dominated by gu — > tbd, giving a Wbbd 
final state after top decay. The resulting multiplic- 
ity distribution looks similar to that of the signal. In 
contrast, the averaged charged particle multiplicity for 
the Wbb and Zbb backgrounds is slightly lower than 
for the signal. This is because the charm quarks from 
stop decay generate on average more charged particles 
than the gauge boson decay products, which contain 
usually only one or three charged particles (the latter 
being due to 3-prong r decays). We have tried sev- 
eral cuts for the charged particle multiplicity. We will 
choose -ZVchargcd > 10 in Sect. IIII El which reduces the 
Wbb and Zbb backgrounds relative to the signal. 

Finally, Fig. [8] shows the number of events as a func- 
tion of the ratio between the pr of the hardest b— jet 
and the p . The signal distribution has a steeper 
fall-off than those for the background processes and 
roughly peaks at one. This is not unexpected, be- 
cause pr(bi) and p T are strongly correlated for the 
signal, due to the relative softness of the second b— jet 
and the c— jets. This means that the fit* pair, whose 
Pt roughly corresponds to the missing pt, essentially 
recoils against the harder 6— jet. From this simplified 
picture we would expect a ratio of < 1. In contrast, 
the single top and ti backgrounds contain more par- 
tons, weakening the correlation between the missing 
Pt and that of any one jet. 



E. Discovery Potential at the LHC 

After we have discussed the basic cuts for our bench- 
mark scenario, we now turn to the discovery potential 
for our bb plus missing energy signature at the LHC. 




5 10 15 20 25 30 35 



pt(bj1)/ptmiss 

FIG. 8: Same as Fig. O but now for the ratio of the pr of 
hardest 6-jet to the missing transverse momentum; more- 
over, we have required that single top and ti events contain 
at least one leptonically decaying W boson. 

We present numerical results for our benchmark sce- 
nario and for heavier stop masses. We show the statis- 
tical significance for the pure QCD prediction of the 
signal and the QCD+EW prediction as a function of 
the stop mass. We shortly discuss the effect of vary- 
ing chargino masses on the statistical significance. We 
also present the statistical significance as a function of 
the stop and neutralino mass. 

From the discussion of kinematical distributions and 
particle multiplicities in the previous subsection, we 
find that the following cuts maximize the significance 
of our signal: 

• -/Vb-jcts > 2, i.e. we require at least two tagged 
6-jets with p T > 20GeV and \rj\ < 2.5. 

• Mepton < 1) i- c - we veto all events with an 
isolated electron or muon with |ry| < 2.5 and 
Pt > 5 GeV. 

• Pr(bi) > 150 GeV, i.e. we require a large trans- 
verse momentum for one of the b— jets. 

• Pt(&2) > 50 GeV, i.e. we impose a significantly 
weaker cut on the second b— jet. 

• p T > 200 GeV, i.e. we demand large missing 
transverse momentum. 

• ^charged > 10, i.e. require at least ten charged 
particles with pt > 2 GeV and \r/\ < 2.5. 

• < 1-6, i.e. we demand that the ratio of 

Ft 

the pt of the most energetic b— jet and p T is not 
large. 

Table IIIII shows the effect of these for our bench- 
mark scenario, assuming an integrated luminosity of 
lOOfb -1 . Numbers are given with a b— tagging effi- 
ciency of 60%. The statistical significance is estimated 
with S = -7=, where S and B are the number of signal 
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events and background events, respectively. We also 
show the ratio S/B. 

We start Table Hill with requiring two tagged b— jets 
with p T > 20 GeV, as well as f T > 200 GeV. With 
these cuts, we already obtain a statistical significance 
of about 49, but S/B is around 0.1. 

After applying the cut on the pr of the hardest 
b— jet, the statistical significance increases to S = 51, 
and the signal to background ratio improves to about 
1/8. Note that this rather stiff cut has relatively lit- 
tle effect on the signal as well as on the Wbb and Zbb 
backgrounds, since the even stiffer f T cut implies a 
rather hard spectrum for one of the b— jets in these 
cases. On the other hand, this cut does reduce the ti 
background significantly. 

After the pr cut on the second 6— jet, the statistical 
significance falls to S = 44. However, for a lower px 
cut on the second 6— jet the tagging efficiency worsens 
[36j , so we keep px > 50 GeV. Note also that this cut 
slightly increases S/B. 

The lepton veto greatly reduces the SM backgrounds 
involving leptonic W decays. Our statistical signifi- 
cance is now S = 66, with a signal to background ratio 
better than 1/4. Most Wbb, ti and single top back- 
ground events passing this cut contain a hadronically 
decaying r lepton, so a hadronic r veto would suppress 
these backgrounds even further. We will come back to 
this point shortly. 

The cut on the charge multiplicity further sup- 
presses the Wbb and Zbb backgrounds, increasing the 
statistical significance to S = 67 and the signal to 
background ratio to nearly 0.3. 

In contrast, the final cut, on the ratio pT(bl)/ft T , 
is quite efficient at suppressing the single top and top 
pair backgrounds. For our benchmark scenario we now 
obtain a statistical significance of S = 69, and a signal 
to background ratio of about 0.35. A good signal to 
background ratio is important for a precise determina- 
tion of the ti — xf~ — b coupling. Furthermore, knowl- 
edge of the systematic error of the SM backgrounds 
at the 10% is sufficient for our benchmark point to 
observe the signal. With an integrated luminosity of 
100 ftT 1 , such a precision is expected to be easily ob- 
tained by the experimental groups [36j . 

As mentioned above, after vetoing events contain- 
ing an isolated electron or muon, most SM background 
events originating from a leptonically decaying W bo- 
son will have a r— jet in the final state. Requiring a r 
veto should reduce the SM background and therefore 
enhance the significance and the signal to background 
ratio. Usually one is interested in identifying, rather 
than vetoing, r jets. The most common method is to 
look for jets above a certain px threshold that contain 
only one or three charged particles, and not too much 
neutral energy; this greatly suppresses backgrounds 
from QCD jets. As an example, assuming a r tagging 
efficiency of 50% for r leptons with px. T > 15 GeV 
and |ry r | < 2.5 we find 7693 background events with 
at least one r in the final state. At the same time, 



the number of signal events stays nearly the same, if 
we assume a small mistagging probability. Therefore, 
the significance is increased to 77, while the signal to 
background ratio increases to 0.43. Vetoing r— jets 
does not reduce the Zbb background; recall, however, 
that this can be subtracted reliably using Z decays 
into electrons or muons. 

The details of the r jet identification algorithm can 
be tuned, allowing to increase the efficiency for cor- 
rectly identifying t— jets at the cost of increasing the 
rate with which QCD jets are misidentified. For ex- 
ample, Ref. [48| finds a 50% r tagging efficiency and a 
~ 98% QCD jet rejection efficiency for privet) ^ 28 
GeV. In the case at hand it would probably be better 
to have a higher efficiency for identifying r— jets, even 
at the cost of more false positives. Recall that at least 
in the absence of ISR and FSR, additional jet activ- 
ity in signal events results from the typically rather 
soft c— quarks produced in t\ decay. The probability 
for misidentifying a c— jet as a t— jet presumably dif- 
fers from that for light quark or gluon jets. Detailed 
knowledge of the detector is thus required for optimiz- 
ing the r identification (or rather, veto) criteria in our 
case. We therefore do not pursue this avenue further. 
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FIG. 9: Significance for QCD contributions only (black 
line) and QCD+EW (red line) as a function of the stop 
mass at the LHC with ^/s = 14TeV. We assume an in- 
tegrated luminosity of 100fb _1 and the mass relations of 
Eq. © and Eq.Q. 

So far, we have discussed the significance for a light 
stop. Now, we turn to larger stop masses. We ap- 
ply the same mass relations for Xi an d Xi as f° r our 
benchmark scenario, i.e. Eqs. © and We use the 
cuts of Table IlIIL In Fig.Hl we present the significance 
for the pure QCD case (black line) and the QCD+EW 
contributions (red line) as a function of the stop mass. 
We assume again an integrated luminosity of lOOfb -1 
and a b— tagging efficiency of 60%. The significance 
for scenarios with the same stop and neutralino mass, 
but with heavier charginos, will lie between these two 
cases; in practice the EW contribution will be ncgligi- 



cut 


Wbb 


single top 


Zbb 


tt 


Signal 


S/B 


s/Vb 


f > 200 GeV, 2 6-jets with 77 < 2.5 
p T {bi) > 20 GeV, p T (6 2 ) > 20 GeV 


6 144 


10 390 


33 440 


179 900 


23 360 


0.098 


49 


p T (bi) > 150 GeV 


5 765 


7824 


27 720 


127 200 


20 760 


0.123 


51 


p T {b 2 ) > 50 GeV 


4 269 


5 476 


19 290 


92 330 


15 360 


0.127 


44 


veto on isolated leptons 


1286 


2 373 


19 290 


32 400 


15 360 


0.278 


66 


# charged hadrons > 10 


1096 


2 227 


15 570 


32 200 


15 020 


0.293 


67 


Pr(bi)/p T < 1.6 


881 


1485 


14 970 


22 030 


13 700 


0.348 


69 



TABLE III: Cut flow for the benchmark scenario of Sect. lIIICl at the LHC with y/s = 14 TeV and an integrated luminosity 
of 100 fb" 1 . Numbers are given for a 6-tagging efficiency of 0.6. 
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FIG. 10: Statistical signal significance at the LHC with 
i/s — 14 TeV as a function of the stop and neutralino 
mass. The red, green and turquoise region corresponds 
to an excess of at least 5a, 3a and 2a, respectively, for 
an integrated luminosity of 100 fb -1 . The chargino mass 
is fixed by Eq. ((4}. The parameter space below the black 
curve is excluded by Tevatron searches [27|, while in the 
region above the blue curve, the stop pair plus single jet 
("mono jet") signal has significance > 5 [3CJ]. The pa- 
rameter region where ti decays into a charm and a neu- 
tralino are expected to dominate is given by the condition 



< m^o + mw + rn b . 



ble (for the given set of cuts) if m~± > 2m; . 

Evidently the EW contribution increases the signif- 
icance by a factor > 2; this is very similar to the 
increase of the total signal cross section found in Ta- 
ble In] After the cuts of Table IIIII the statistical 5tr 
discovery reach extends to 260 GeV (210 GeV) with 
(without) the EW contributions. 

In Figure [TOJ we show the significance in the stop- 
ncutralino mass plane for an integrated luminosity of 
100 fb -1 at s/s = 14 TeV. We assume the mass relation 
of Eq. (HJ). The parameter space below the black curve 
is excluded by Tevatron searches at the 95% confidence 



level [27] , whereas the parameter points above the blue 
curve allow detection of a monojet signal with at least 
5cr significance [30l ]. In our analysis, we assume that 
the light stop decays dominantly into a charm quark 
and the lightest neutralino. This will generally not be 
the the case if the mass relations m; > m c o + m r and 
mj-j < m^o + mw + m& do not hold. 

We first observe in Figure [TO] that the significance 
increases with decreasing stop-ncutralino mass differ- 
ence [HD]- Recall that we require pribi) > 150 GeV. 
This implies that the t\t\ pair usually has a rather 
large transverse momentum, which in turn leads to 
significant boosts of the stop squarks. As a result the 
charm jets tend to go into a direction close to that 
of the stop pair; i.e. the charm jets tend to reduce 
the missing px- Scenarios with softer charm jets, i.e. 
smaller t\ — Xi mass splitting, therefore pass thef T cut 
more efficiently. A similar argument also holds for the 
significance of the monojet signal, where the missing 
energy is is required to have fax > 1 TeV; moreover, a 
veto against additional jets is imposed [30j . 

Our results suggest that a discovery of stops via 
stop pair production in association with two b— jets is 
possible as long as < 270 GeV. In the case of non- 
observation of any signal one would be able to exclude 
the parameter space at 2a for stop masses of up to 
340 GeV; cf. the turquoise region in Fig. [TOJ Our sig- 
nal is also visible in regions of parameter space where 
the monojet signature produces no significant excess 
over the SM backgrounds [54| ■ This is partly because 
our process is additionally enhanced by the EW con- 
tributions, which do not enter significantly the mono- 
jet production process. We note however, that differ- 
ent cuts for the monojet signal than those of Ref. [30| 
might increase the monojet discovery reach [55j. 



IV. POSSIBILITIES TO FURTHER 
OPTIMIZE THE CUTS 



So far we have only applied cuts that can be used for 
all combinations of stop and LSP masses, i.e. we have 
not tried to optimize the cuts as function of or 
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to^o. This simplifies the statistical analysis: if we try 
several combinations of cuts, the chance of finding an 
upward fluctuation of the background increases. On 
the other hand, additional cuts could help to increase 
not only the statistical significance, but also the pu- 
rity of the sample, i.e. the signal to background ra- 
tio. Note that Table Hill implies a signal to noise ratio 
of about 0.025 at the limit where the signal becomes 
statistically significant, i.e. for S/yB = 5. If we re- 
quire S/B > 0.1, the search reach would decrease from 
about 270 GeV to 190 GeV. After imposing the r veto, 
this would improve again to about 200 GeV, which is 
however still much smaller than the statistical reach. 

Recall that our background is dominated by Zbb and 
ti events. The systematic uncertainty on the former 
is very small, since one can analyze events with iden- 
tical kinematics where the Z decays into an e + e~ or 
pair. Of course, ti events will also be analyzed 
in great detail by the LHC experiments; nevertheless 
some extrapolation into the signal region will presum- 
ably be required to determine this background, which 
will introduce some systematic uncertainty. 

Note also that the statistical 5er discovery limit after 
the cuts of Table IIIII still corresponds to about 1 000 
signal events. Refined cuts could therefore also help 
to further increase the mass reach. 

In this section we therefore first discuss some cuts 
that might help to suppress the ti background for sce- 
narios with small t\ — \\ mass splitting. We then 
discuss possibilities to identify the c— jets in the sig- 
nal, which can be used to suppress all backgrounds, 
and also to distinguish our scenario from other SUSY 
processes with a similar final state. 



A. Further cuts for small mass splitting 

If mjj — m^o is small, the c— jets will usually be too 
soft to be detected. In contrast, many ti background 
events will contain a hadronically decaying W boson 
in addition to the two required 6— jets. We saw in 
Fig. [7] that this leads to a higher charged multiplicity 
of these background events. One can cut on several 
(related) quantities to suppress this background (and, 
to a lesser extent, the single top background): 

• A jet veto could be imposed. Recall, however, 
that the signal is almost exclusively due to gluon 
fusion, which implies strong initial state radia- 
tion. In order not to unduly suppress the sig- 
nal it might therefore be preferable to only veto 
events that contain two (or more) additional 
hard jets. 

• An upper limit on the total transverse energy 
could be imposed; probably the two b— jets 
should not be included here, i.e. one should only 
sum over particles (or calorimeter cells) that are 
not part of these two jets. 



• An upper limit on the charged particle multiplic- 
ity could be imposed. For example, requiring 
-^charged < 20 would reduce the ti background 
by ~ 30% with little loss of signal. 

B. Charm tagging 

Observing two soft charm jets in addition to two 
hard b— jets and large missing energy would give strong 
evidence for our scenario. However, soft charm jets 
are probably difficult to identify at the LHC. The un- 
derlying event, QCD radiation and possibly overlap- 
ping events (which we did not include in our simula- 
tion) generate substantial hadronic activity in our sig- 
nal events and may overshadow potential charm-jets. 
That said, we can still expect some reasonably ener- 
getic c— jets for the light stop scenarios if the i\ — Xi 
mass splitting is not too small. For example, for our 
benchmark point and after applying all the cuts of 
Table IIIII we expect about 45% of all signal events to 
contain at least one c— jet with px > 50 GeV [56| . 

Still, we have to discriminate the charm jets from 
light-flavored and gluon jets. In Ref. [3(|, the authors 
use the jet mass and charged particle multiplicity in- 
side the jet to discriminate between c— jets and light 
flavor or gluon jets. They find c— tagging efficiencies 
of > 50% for mass splitting > 10 GeV, but with a 25% 
mis-tagging probability. In the case at hand it might 
be better to employ a more sophisticated algorithm, 
e.g. based on a neural net, that attempts to decide 
whether the hadronic activity not associated with the 
two 6— jets is consistent with containing two c— jets, 
rather than identifying the c— jets in isolation. 

Identifying at least one c— jet would suppress the 
Wbb, single top and Zbb backgrounds. Moreover, it 
would distinguish the stop pair production we are con- 
sidering from the pair production of light sbottom 
squarks (b\), followed by b\ — > bxl decays, which also 
leads to a final state with two b— jets and large missing 
Pt, but without c— jets. However, half of all ti events 
containing a hadronically decaying W boson also con- 
tain a c— jet (from a W + — > cs decay or its charge 
conjugate); c— tagging will therefore probably not be 
as efficient in reducing the ti background. 

V. SUMMARY AND CONCLUSION 

In this article, we analyzed stop pair production in 
association with two b— quarks in the co-annihilation 
region, where the t\ — Xi mass difference is small, 
so that simple stop pair production does not yield 
an observable signal. Assuming a light higgsino-like 
chargino beside a stop NLSP and neutralino LSP, the 
leading order QCD as well as mixed QCD-EW contri- 
butions are taken into account. We discussed in some 
detail the dominant diagrams for our signal process. 
We simulated the signal process and dominant SM 
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background processes. For a benchmark scenario with 
mi = 120 GeV, m,« = 100 GeV and m~± = 140 GeV, 
we described all important kinematic distributions and 
developed a set of cuts that help to isolate our signal 
process. The most effective cut is a veto on isolated 
leptons. We showed that we can have a significance 
of S = 70 for our benchmark point, assuming an inte- 
grated luminosity of 100fb _1 at ^/s = 14TcV; alter- 
natively, a signal with 5cr statistical significance can be 
obtained with just 0.5 fb _1 of data. For fixed xf — ti 
and t\ — Xi mass differences, the same set of cuts allows 
to detect the signal at > 5a statistical significance for 
mjj < 270 GeV. Note that the significance of our sig- 
nal actually increases when the f i — X? mass difference 
is reduced, making it complementary to searches for 
direct stop pair production. The inclusion of the EW 
diagrams greatly enhances the total cross section and 
more than doubles the statistical significance. 

However, we saw that for > 190 GeV our cuts 
leave us with a relatively poor signal to background 
ratio. This can be improved by vetoing r— jets, i.e. 
hadronically decaying r leptons. We also briefly dis- 
cussed additional cuts on the hadronic activity not 
associated with the two hard b— jets that could help 
to further suppress the ti background, which proba- 
bly will have larger systematic uncertainties than the 
Zbb background. These cuts should be tailored to the 
t\ — Xi mass difference, and perhaps also to the over- 
all mass scale (which affects the amount of QCD ra- 
diation); we have therefore not attempted to inves- 
tigate this systematically. We also briefly discussed 
the possibility to use charm tagging to further sup- 
press the backgrounds, and to distinguish our process 



from other SUSY reactions yielding two b— jets and 
missing px- Charm tagging, and r vetoing, depend 
quite strongly on details of the detector performance; 
a quantitative analysis of their efficiency is therefore 
best left to our experimental colleagues. 

One motivation for analyzing this final state is that 
the mixed QCD-EW production channels might al- 
low to probe the t\ — Xi — b coupling, thereby allow- 
ing for the first time to check a SUSY coupling rela- 
tion involving Yukawa couplings; we intend to inves- 
tigate the feasibility of determining this coupling in a 
later publication. This will require that the t\ mass is 
known. Determining both the mass and the EW cou- 
pling should be easier if the production of stop pairs 
recoiling against a photon or a very hard jet [3 Of also 
leads to an independent statistically significant signal, 
which does not depend on the t\ — x^ — b coupling. 
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We do not consider additional CP-violating phases 
(compared to the SM one) in this paper. They are 
unimportant for our collider process. 
We have calculated the Xi relic abundance for various 
of our scenarios with the help of micr0megas2.4 [2p( ]. 
The corresponding splitting function is proportional 
to x 2 + (1 — x) 2 , i.e. it peaks at x = and x = 1. 
Figs. [7] and [H] show only single top and ti events with 
at least one isolated electron or muon in the final state 
or events with at least one hadronically decaying tau. 
We employed the lepton cuts of Sect. IIII El 
Small deviations from this behaviour are due to sta- 
tistical fluctuations of the signal. 

We note however that in contrast to our work, Ref. [3(| 
included systematic uncertainties equal to \j7B; see 
Ref. [3(|] for more details. The monojet discovery reach 
will therefore increase if one uses solely S/VB as an 
estimator for the significance. 

We did some estimates for the monojet signal and the 
most important backgrounds and found that less hard 
cuts on the monojet pr and on the far lead in general 
to larger significances. We plan to present our results 
in a future publication. 

Recall that we use R = 1 in our jet definition. This 
allows a significant contribution from the underlying 
event and from unrelated QCD radiation to the jet 
energy. In the absence of these effects, only ~ 16% of 
all signal events contain a c— jet with pr > 50 GeV. 



